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It is shown that the nonequihbrium steady-state of a phase-coherent conductor can be described 
by a generalized thermodynamic potential based on the concept of partial densities of states. This 
is possible due to the fact that dissipation takes place only in the contacts for the emitted carriers. 
O ' Long-range Coulomb interaction is included, and charge conservation and gauge invariance are 

2^ \ satisfied. The theory is illustrated for a mesoscopic capacitor with leakage. 
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fH ' The density of states (DOS) is the basic quantity from which equilibrium thermodynamic properties of a system 
I , are derived. Once the DOS is known, one can construct an appropriate thermodynamic potential which has to be 
^ • minimized in order to find the equilibrium state. Due to a production of entropy in driven systems, on the other 
hand, nonequilibrium states can in general not be obtained by minimizing a thermodynamic potential. In this work 
I I show, however, that generalized thermodynamic potentials still exist in mesoscopic conductors |^ in the presence 
of transport, provided that dissipation takes place only in the reservoirs In order to treat the nonequilibrium 
case, the concept of the DOS is generalized by introducing partial densities of states |^ which contain the information 
from which reservoir the particles are injected. According to Biittiker the specific partial DOS used below are 
'"^ , called injectivities. The injectivities play an important role in the theory of time dependent and nonlinear electrical 
P5 transport in mesoscopic conductors One can show that the injectivities are related to the dwell-time of the 

O ' particles in the conductor The injectivities appear naturally within the scattering approach to conduction and 

I ^ 1 1 can be expressed in terms of a scattering matrix Q . Equilibrium thermodynamics in terms of scattering matrices has 
been formulated in a different context by Dashen, Ma, and Burnstein Avishai and Band ||^ discussed the relation 
T-H , between the scattering matrix and the total DOS of a one-dimensional system. 

^ ■ This paper is organized as follows. First, the definition of the injectivities is recalled for a system of non-interacting 
particles. Secondly, the generalized thermodynamic potential including the Coulomb interaction is constructed for a 
' phase-coherent multi-terminal conductor. The charge distribution is then determined by minimization of this poten- 
] tial. Finally, the formalism is applied to a symmetric capacitor with tunneling between the capacitor plates and it is 
• shown that the result is in accordance with Ref. . 
\^ ' Consider a phase-coherent conductor connected via contacts to reservoirs a = 1, TV of non-interacting electrons. 
On ■ Such a conductor is characterized by a unitary single-particle scattering matrix which can be decomposed into sub- 
matrices Sap- The indices a and /? label contacts, and the dimensions of the submatrices are equal to the numbers of 
channels in the associated contacts. The first and the second index of Sap correspond to out-going and to incoming 
y ' particles, respectively. The scattering matrix is a function of the energy E of the scattered particle, and is a functional 
' , of the single-particle potential eU{x), where e is the electron charge and where U{x) is the electric potential. Denoting 
^ the trace of a matrix by 'Tr', one can write the transmission probability of a carrier from contact /3 to contact a as 
O Taf3 = Tr(s|j^s„^). The reflection probability at contact a is Ra = Tr{sj^^Saa)- At thermodynamic equilibrium, the 
O ■ response of the particle density to a variation of the Fermi energy is characterized by the local DOS which can be 
^ I expressed in terms of the scattering matrix elements and its functional derivatives with respect to the potential |^J^ 
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which are the partial DOS associated with particles injected at contact (3. In a nonequilibrium situation where 
the electrochemical potential in a single reservoir is changed, the response of the particle density is given by 
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5n{x) = {dn{x, f3)/dE)dfj,p. A decomposition of the local DOS into injectivities leads to the following picture. The 
total sample is decomposed into N subsystems. Subsystem f3 consists of all those scattering states which are associated 
with particles injected at contact (3. The local DOS of subsystem f3 is given by the injectivity dn{x,(3)/dE. In a 
nonequilibrium steady-state, these states are filled as if subsystem (3 were in equilibrium with reservoir (3. Dissipation 
takes place in the reservoirs only for the out-going particles which are thermalized to the Fermi distribution of this 
reservoir. As long as scattering in the conductor is elastic, dissipation does not affect the states inside the sample. 

For simplicity, I assume zero temperature and consider a discretized version of a mesoscopic conductor. Importantly, 
all nearby conductors and gates are included in this model. The whole system consists thus of regions flk {k = 1, M) 
with electrostatic potentials C/fe, charges qu, and injectivities Dka- For later convenience, the injectivities are written 
in the form 

Dk^ = [ Sx , (3) 

which have the dimension of a capacitance. The local DOS of a single region Q.k and the total DOS of the sample are 
Dk = ^f^a ^'^'^ ^ — J2ka ^ka, respectively. If eVa'^'^^ denotes the equilibrium electrochemical potential of contact 
a, the part of the charge on conductor k which is injected from contact a for a variation AVa — Va — vlf'^'^ of the 
voltage becomes 

fta = / ° dVDuo.{V) . (4) 

Jvi""'' 

The total charge in fi^ is qk = X^a^fc"- The energy E^ of the noninteracting system is given by the sum of the 
energies over all single-particle states 



ED=y2 dVVDUV) , (5) 



where the energy scale is defined such that the equilibrium system has zero energy, E^^"^^ ~ 0. After an expansion 
of Eqs. (^) and (||) to second order in AVq, one can eliminate the AVa and express the energy as a function of the 
partial charges qka 



In order to include the long-range Coulomb interaction, it is convenient to introduce a geometric capacitance matrix 
Cki for the regions flk, which is determined by the Poisson equation. An arbitrary charge distribution {qi} induces 
electrostatic potential shifts Uk given by 
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Uk=J2^k'li + Uo . (7) 
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Note that a global voltage shift Uq is always a solution of the Poisson equation. The Coulomb energy of the charge 
distribution is thus 

Ec{{qka})^^^qkaC^l^qil3 + Uo^qka ■ (8) 
kaLj3 ka 

At zero temperature, the free energy is equal to the total energy of the closed system and is given by the sum of the 
kinetic energy Ed and the Coulomb energy Ec 



EiUka}) = ^(t/o + V^^''^)qka + If^-C^aiplip , (9) 

where the following matrix is introduced 
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Ckai0 = D^aSkiSaf^ + ■ (10) 
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Since the open system is appropriately described by the grand-canonical ensemble, the potential which must be 
minimized is 

E{{qkc.}) - J2 ^"^'^^ = ■ (11) 

ka. 

Variation with respect to the charges qka yields finally a set of Af x equations 

E^fc"W = ^^--t^o , (12) 

1/3 

where AVa — Va— vif^^ corresponds to the (electrochemical) voltage shift in contact a. The still free global shift ?7o 
of the electric potential is determined by the additional condition of charge conservation, 



^gfca =0 . (13) 
ka 

Thus, C/q can be interpreted as a Lagrange parameter associated with the condition (|l3|). Charge conservation is a 
general property of a complete set of conductors connected to electron reservoirs, since electric fields are fully screened 
in the reservoirs j^. In order to solve Eq. ( |l^ ) for the qkai the matrix Ckaip is introduced as the inverse matrix 
of C^^ip. Note that this quadratic matrix acts on M x A^-dimensional vectors q^a (see also the example below). 
Combination of Eqs. (|l2|) and (|l^) yields 

° - 7< ■ (1^) 
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Once the partial charges qka are known, one can calculate the electrochemical capacitance which relates the static 
charge distribution q^ — qka to the voltage shifts AVa in the contacts, 

_ dqi _ J2k/3mjnsi^iPkaCmjnS — Clf^nsCm-ika) 

The nonequilibrium electric potential follows from the charge distribution with the help of Eq. (^ . 

An important consequence of charge conservation is that the electrochemical capacitance matrix satisfies C^./q = 
1^. Additionally, the result is gauge invariant, J2a ^t^M = [||, i.e. it does not depend on a global voltage shift. 
The global voltage shift is absorbed by the constant Ua and does not change the total energy of the system since the 
total charge vanishes. 

As a simple example, consider a symmetric two-terminal capacitor with tunneling between the capacitor plates 
(leakage). Assume a single one-dimensional open channel with a transmission probability T ~ \ — R. Equivalently, 
this system describes a symmetric quantum point contact with a single open channel, or a one-dimensional conductor 
containing a symmetric impurity. This system has already been discussed in Ref. Q in the context of time-dependent 
transport. I show that the thermodynamic treatment introduced in this work is in accordance with the results of Ref. 
1^. Within a semiclassical approximation, the injectivity from reservoir 2 to plate 1 is proportional to one half of the 
DOS, Di {— D2), of plate 1 and to the transmission probability, T. The factor one half occurs since only particles with 
a velocity in direction to contact 1 contribute to this partial DOS. Due to symmetry, it holds D12 — D21 — DiT/2 
which imphes Dn = D22 = Di - D12 = Di{l - T/2) The diagonal elements A = C^^ = C22 and the off- 

diagonal elements B = C^^ — ^21 of the inverse geometrical capacitance matrix are given by A = C\\I[C\]^ — C12) 
and B — —Ci2/{Cii — 6*^2)- The effective (charge conserving) geometric capacitance Co between the two capacitor 
plates is given by Co — (Cn — Ci2)/2. In order to represent the matrix Ckaip in a simple way, it is convenient to 
abbreviate Ki = D^^ + A and K2 = + ^- The equations (n2) read then 



Ki A B B \ / qii\ / AVi-Ui 



B B K2 
A K2 B 








AVi - Uo 
AV2 - Uo 



(16) 



B B A KiJ \q22J \AV2~U0 



By taking charge conservation ( |l^ ) into account, these equations can be solved for the q^a- From the charge qi 
^92 = 911 + qi2 on the capacitor plates one finds the electrochemical capacitance 
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dqi dqi R 



This result states that, in the case of vanishing transmission, the electrochemical capacitance is the geometrical 
capacitance in series with quantum corrections given by the DOS of the plates jj]. In the macroscopic limit where the 
DOS diverges {D — > oo), one recovers the pure geometric capacitance, = Cq. On the other hand, for increasing 
transmission the capacitance decreases and vanishes for T = 1. Clearly, without reflection there is no charge dipole 
at all. A more extensive discussion is provided by Refs. It should be clear that transport properties such as, 

e.g., the current cannot be calculated directly from such a generalized thermodynamic potential which contains only 
the information on the charge (or potential) distribution. 

In conclusion, I have shown that it is possible to describe the stationary nonequilibrium charge distribution of 
a mesoscopic phase-coherent conductor with a generalized thermodynamic potential. The theory is based on the 
concept of partial densities of states (injectivities). The result includes the long-range Coulomb-interaction and is 
charge conserving and gauge invariant. To illustrate the theory, I used the simple approximation of a descretized 
version of the mesoscopic sample. However, in principle the generalization to a field theory, e.g., in the framework of a 
density functional theory, is straightforward. The generalization to finite temperatures and to other thermodynamic 
questions (e.g., mechanical properties of the system), remains an interesting future task but does not require new 
concepts. Furthermore, I considered only the deviation from the equilibrium state described by equal electrochemical 
potentials. But this is not a restriction since the equilibrium reference state can always be constructed from a 
thermodynamic potential. Although transport properties cannot be extracted directly from the generalized potential, 
it must be emphazised that the knowledge of the nonequilibrium state is crucial for the determination of the currents 
beyond linear dc-response [^|-^| JlC|Jl^ . 
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